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ABSTRACT 

^ Aims. The ASTEP (Antarctica Search for Transiting ExoPlanets) program was originally aimed at probing the quality of the Dome C, 

Antarctica for the discovery and characterization of exoplanets by photometry. In the first year of operation of the 40 cm ASTEP 400 
telescope (austral winter 2010), we targeted the known transiting planet WASP- 19b in order to try to detect its secondary transits in 
^ the visible. This is made possible by the excellent sub-millimagnitude precision of the binned data. 

I Methods. The WASP- 19 system was observed during 24 nights in May 2010. Once brought back from Antarctica, the data were 

O processed using various methods, the best results being obtained with an implementation of the Optimal Image Subtraction (OIS) 

^ algorithm. 

(j^ Results. The photometric variability level due to starspots is about 1.8% (peak-to-peak), in line with the SuperWASP data from 2007 

Ci (1.4%) and larger than in 2008 (0.07%). We find a rotation period of WASP-19 of 10.7 + 0.5 days, in agreement with the SuperWASP 

determination of 10.5 + 0.2 days. Theoretical models show that this can only be explained if tidal dissipation in the star is weak, 
i.e. the tidal dissipation factor 2* > 3 x 10^. Separately, we find evidence for a secondary eclipse of depth 390 + 190 ppm with 
a 2.0cr significance, a phase consistent with a circular orbit and a 3% false positive probability. Given the wavelength range of the 
observations (420 to 950 nm), the secondary transit depth translates into a day side brightness temperature of 2690^220 ^^^^ ^^^^ 
measurements in the and K bands. The day side emission observed in the visible could be due either to thermal emission of an 
extremely hot day side with very little redistribution of heat to the night side, or to direct reflection of stellar light with a maximum 
geometrical albedo Ag = 0.27 + 0.13. We also report a low-frequency oscillation well in phase at the planet orbital period, but with a 
lower-limit amplitude that could not be attributed to the planet phase alone, and possibly contaminated with residual lightcurve trends. 
^T) Conclusions. This first evidence of a secondary eclipse in the visible from the ground demonstrates the high potential of Dome C, 

Antarctica for continuous photometric observations of stars with exoplanets. These continuous observations are required to understand 
^T) star-planet interactions and the dynamical properties of exoplanetary atmospheres. 

\ I Key words. Instrumentation: photometers - stars: individual: WASP-19 - (stars:) planetary systems - (stars:) starspots - Planets 
^ and satellites: atmospheres 

^ 1 ■ Introduction ever, in spite of a very rapid increase of the number of exoplan- 

, . r . , , , . ets for which these detections have become possible, our under- 

The observation of primary and secondary planetary transits at ^^^^^j ^^^^^ atmospheres remains sketchy at most. This is 
multiple wavelengths is a powerful tool to characterize the at- ^^^^^^ ^^^^ detections have been done at infrared 

mospheres of exoplanets, understand their energy budget, atmo- ^^^eie ths and have therefore concerned the thermal part of 
sphenc dynamics and thermal evolution. The first detections of spectrum, rather than also the visible part which is 

a secondary transit fDemingetal.'2005V of phase variations of a ^^^^^^^j^j ^^^j^ ^^^^^^ irradiation absorbed by these 
non-transiting planet (Harrington et al. 2006) and of a transiting atmospheres. The situation has been changi ng thanks to obser - 
planet (Knutson et al. 2007), confirmed theoretical predictions ^^^-^^^ ^ th e space telescopes MOST (e.g. | Rowe et al. | 2008 |>. 



Of a Strong day-night temperature contrast and power^^^ zonal ^^^^^ ^ ' Alonso et al. |[2 009aF, " Snellen e t aL || 2QQ9 tW 
wmds m close-m exoplanets ([Showman & Guillot|2002|). How- ^^^^^^ (e. g. |Bomc]d et al. 2009 , Kipping & Lkos |[2m]i: 



Send offprint requests to: L. Abe, but progress has been slow due to their intrinsic limitations in 

e-mail: Lyu.Abe@unice.fr accessible targets. 
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We present near-continuous observations of the WASP- 19 
star-planet system acquired at visible wavelengths from the 
ASTEP 400 telescope at the Concordia station, Dome C, Antarc- 
tica. As described by |Crouzet et al.| (|2010 ), the DomeC has 
excellent meteorological conditions and near-continuous obser- 
vations are possible for approximately three months centered 
on the Southern winter solstice (June 21st). This makes it an 
ideal site for long-duration observations such as those required 
to characterize the phase curves of exoplanets. 

The target, WASP- 19, is hosting one of the exoplanets with 
the shortest orbital period known thus far. It was discovered 



as part of the WA SP survey in the southern hemisphere ( Hebb 
eraL||2QTQ| [HTQ| hereafter). Its visible magnitude of 12.3 and 
declination of about -45.66° made it a perfect target for ob- 
servation from Concordia. The atmospheric properties of exo- 
planet WASP- 19b are beginning to be well-characterized in the 
infrared thanks to observations of its secondary eclipses at 1.62 
and 2.09 ^im us ing VLT/HAWK-I (Anderson et al. 20Tq1 



|son et al.||2010[ respectiv ely), at 3.6, 4.5, 5.8 and 8.0 |Lim using 
Spitzer (l Anderson et al.||2011| ), and recently in the z' band at 
0.9 ^tm with the NTT/ULTRACAM ( [Burton et al. ] [20T2l ). The 
interest in the system is further enhanced due to its short orbital 
period (0.79 days) and Jupiter-like mass (1.17 Mjup) which im- 
plies that its survival is critically dependent on the efficiency of 
tidal dissipation in its solar-type host star (0.97 M©). Any con- 
straint on the orbital properties of the system, in particular or- 
bital eccentricity and spin period of the star can help us to better 
understand its dynamical evolution, fate, with potential implica- 
tions for the entire population of close-in exoplanets. The article 
is organized as follows: We first present the observational setup 
and processing of the data. The results are presented in Sect.|3] 
In Sect.|4j we discuss how these results compare with other well 
characterized exoplanetary atmospheres, and how these mea- 
surements can be applied to understand the dynamical evolution 
of the systems. 



Table 1. The ASTEP 400 WASP- 19 observations parameters and statis- 
tics. 



Data 


Value 


Notes 


Average temperature 


-68±7°C 




Exposure time 


130s 


single frame 


Frame periodicity 


150 s 


median 


PSE EWHM 


4 arcsec (4.3 pix) 


median 


Max. sun elevation 


-10° to -6° 


daily 


Instrument Duty cycle* 


61% 




Overall data duty cycle 


58% 


~14hr daily 


Science only data duty cycle 


48% 





* fraction of time of instrument operation including science and 
calibration data, telescope and systems overheads. 
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Fig. 1. Transmission in the ASTEP 400 camera as a function of wave- 
length (in nm). The vertical lines indicate the range of wavelengths 
within which 65% of the flux of a T^^ - 5500 K blackbody is transmit- 
ted, i.e. 575 to 760 nm. 



2. Observations and data processing 

The ASTEP (Antarctica Search for Transiting Ex o-Planet) pro- 
gram comprises the ASTEP-South 10 cm refractor ([Crouzet et al.| 
[2010^ and the ASTEP 400 40 cm telescope ( Daban et al.||2010| ) 
that was used for the WASP- 19 observations reported in the 
present paper. The goal of the observations was not only to try 
to detect a signature of the secondary transit, but also the possi- 
ble reflected light phase variations of the planet, and the stellar 
luminosity variation reported in HIO^ We therefore observed the 
target star continuously for 24 days. 



2.1. Observation log 

WASP- 19 was one of our first ASTEP 400 program targets. It 
was observed between April 30* and May 23^^ 2010 (UTC). 
We summarize the basic observation parameters and statistical 
information in Table[T] The information concerning duty cycle 
(instrument and data recording) is measured relative to the over- 
all time span of the WASP- 19 campaign. 

The temperature ranged from -78°C to -44° C with an average 
value of -67° C. On the instrument side, we already noticed ice 
forming on the primary and sometimes on the secondary mirror. 
At that time we did not have any permanent defrosting system 
installed, so a daily manual defrosting operation was mandatory. 
Icing on the optics did of course aff'ect the absolute photometry, 
but since it aff'ected all PSFs in the field, we could obtain good 
lightcurves by careful data processing (see below). 



2.2. ASTEP 400 transmission 

The ASTEP 400 telescope consists of a 400 mm Newton tele- 
scope designed to withstand the harsh conditions in Antarctica 
and a sophisticated field corrector. The beam is reflected by 
2 aluminum-coated mirrors, before being sent to the thermal- 
ized field corrector where it first crosses a double lens designed 
to maximize both thermal isolation and image uniformity. The 
beam is then split by a dichroic so that the blue part is used for 
guiding and the red part of the light is used for the science cam- 
era. On its way to the science camera, the beam passes again 
through three lenses plus the CCD glass window. The science 
camera is an ELI ProHne with a KAE 16801E, 4096 x 4096 pix- 
els front-illuminated CCD and a 16 bits analog to digital con- 
verter (see [Daban et"aL]|2010| for a complete description of the 
instrument). 

Eigure [T] shows the final transmittance for the science im- 
age as a function of wavelength calculated using Zemax and 
technical data obtained from the constructors for each element 
and the CCD camera. Most of the signal is transmitted in the 
red. In order to estimate more precisely our wavelength cover- 
age for science images, we assume that transmission stops be- 
low 300 nm and above 1000 nm. We then calculate the range of 
wavelengths that correspond to a 65% and a 95% transmission 
of a Tqq = 5500 K blackbody, corresponding to that of WASP- 
19. The result is 575 to 760 nm and 420 to 950 nm, respectively. 
(In the case of a uniform irradiation, we find that 65% of the flux 
would be transmitted between 580 and 780 nm.) 
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2.3. Data reduction 

Each of the ASTEP 400 frames is pre-processed using standard 
bias and dark corrections. Elat field correction turned to be 
a little more complicated sin ce we had to deal wi th the sky- 
concentration eff'ect (see e.g. [Andersen et al.|[T995] ). We used 
satellite trails and star brightness statistics to infer the eff'ect of 
the sky-concentration, and confirm that the bell-shaped halo at 
the center of our images was a purely additive component. We 
removed it by performing a 32 x 32 median smoothing on a 
256 X 256 reduced copy of the frame. The enlarged version of 
this smoothed image is then subtracted from the original frame. 
This process had to be applied to each individual frame since 
the sky-concentration changed slightly over time. The flat field 
frames were obtained with the same process. 

The ASTEP 400 data pipeline exists under two difl'erent im- 
plementations: a custom IDL code using classical aperture pho- 
tometry routines from the well known IDL astronomical library, 
and the Miller-Buie implementation of the Optimal Image Sub- 
traction algorithm (QIS, [Miller et al.|200 8 ). The first implemen- 
tation is used for processing the data on-site at the Concordia 
station, whilst the second one was developed later. 

The OIS core algorithm (see Appendix[A| is similar to the 
ISIS algorithm ( [Alard & Lupton||1998| ) with the diff'erence that 
in ISIS the convolution kernels are made of combinations of 2- 
D Gaussian functions, whilst in the OIS, these kernel functions 
are completely arbitrary (the optimal ones with respect to the 
subtraction residual). The superiority of IS IS -type algorithms 
over aperture photometry has been clearly demonstrated for ex- 
ample in |Montalto et al.| ( |2007[ ). The OIS requires much more 
computational time, but also significantly improves the scatter 
in all our lightcurves compared to aperture photometry by a fac- 
tor -1.5 in RMS see Appendix ] A. 5 [ There is also an improve- 
ment of the noise structures (red noise) that makes the OIS data 
globally cleaner. Specifically for the WASP- 19 target which is 
amongst the brightest stars in that field, the gain is 20%. Quite 
surprisingly very few papers (e.g. |Gillon et"ar||2007| ) report on 
the systematic photometric biases that is intrinsic to the image 
subtraction methods. A description of this eff'ect and a way to 
compensate is discussed in Appendix[A[ In the case of WASP- 
19, the RMS gain compared to aperture photometry is about 20% 
(see folded lightcurve comparison in Appendix ] A 5 1 ). 

The OIS is intrinsically removing all global variations for all 
stars. The remaining systematics are calibrated using the follow- 
ing methodology: about 3000 stars are pre-selected as not be- 
ing noticeably varying (using an unsharp mask filter on the raw 
lightcurves and keeping stars with less RMS scatter). Amongst 
these 3000 remaining stars, only 10 stars are kept as references 
for each of the analyzed star. The selection criteria consists in 
keeping those 10 stars that minimize the most the RMS scatter 
after the lightcurve is normalized. We have compared our cali- 
bration procedure with other methods such as SysRem ( [Tamuz, 



Mazeh, & Zucker]|200 5), or the method described in ( [Mazeh 
et al. 2009| ) without noticeable diff'erences. 



Prior to applying primary transit and secondary eclipse 
model fitting (see Sect.|3.2.l] and |3.3|), the obtained lightcurve 



(Eig.[2]-top) needed to be filtered to remove the stellar varia- 
tions. A calibration lightcurve (excluding primary transit win- 
dows) was obtained by a temporal median filter, that was then 
smoothed by applying a Eourier- space filter. This smoothed 
lightcurve was subtracted from the original one. The obtained 
residual lightcurve (flattened) was median-filtered again. The 
temporal width of this last filter was obtained by running a simu- 
lation where the primary transit was simply replaced by artificial 



500ppm-deep transit events. The optimal filter width was found 
to be 6.25 hr, as it was selected to maximize the number of arti- 
ficial eclipse detections. The final calibration lightcurve is com- 
pleted by interpolating the out-of-transit smoothed lightcurve 
at the primary transit locations. The "flattened" lightcurve 
(Eig.[2]-bottom) is obtained by subtracting the smoothed vari- 
ability lightcurve to the original data. 

2.4. Data quality 

Eor each night, the pink noise - gaussian and corr elated noise 
combined - was analyzed with a variant of the |Pont et aL 



( ,2006) prescription, taking into account the preprocessing of the 
lightcurves. As described in ( [Gillon et al.||2006| ), it consists in 
computing the variance of the pink (white + red) noise on the 
lightcurves, and on a smoothed version of them (the bin size 
corresponding to the duration of the transit). Each noise vari- 
ance follows the relation cTp^j^j^ = ^wMte 

/N + (T^g^ where N is 
the number of points used for smoothing the data. The solution 
of this system of two equations are the white and red standard 
deviations (Xwhite and (Xwhite- The difl'erent nights were found 
to vary in stability, with the hour-timescale pink noise at best 
430 ± 65 ppm, median 700 ppm and worst 955 ± 250 ppm, com- 
parable to very good ground-based transit photometry. The un- 
certainties on these values were evaluated by generating artificial 
data with corresponding white and red noises statistics. The red 
noise uncertainty is 200 ppm in average. 

The result of this analysis (Eig.|5]) shows that 3 nights are 
"not-so-good", with a pink noise over 1 hour around 0.1 % (night 
7 has too few points to compute relevant statistics, but was kept 
as it adds information for the stellar variability). Most nights 
are "good", with a pink noise over 1 hour around 0.7mmag, 
and a few "very good" nights, with a pink noise over 1 hour 
around 400 ppm. These numbers are typical to current ground- 
based photometry. The best ground-based transit photometry 
now reaches 0.4 mmag (KeplerCam, VLT infrared), the "good" 
nights are also around 0.7, and more casual observations are 
above 1 mmag. 

During the first night, we noticed that there was a switch of 
guiding star from the good one to another one. Taking into ac- 
count this first night produced odd results (strong trends) when 
processing for the photometry with the OIS pipeline. In addition, 
the fraction of observing time when the telescope was tracking 
the correct guide star was too short to aff'ect our result signif- 
icantly so we deliberately discarded all data points from that 
night. In the rest of this paper, all results are obtained from the 
23 remaining nights. 



3. Results 

3.1. Stellar variability 

Before the ASTEP 400 campaign, the only observations with a 
lo ng-en ough duration to extract the stellar spin period were those 
of |H10| who found a periodicity of 10.5 + 0.2 days. Our data also 
exhibit a periodic oscillation with an amplitude of about 0.9% as 
well as a period of 10.7+0.5 days, compatible with the 2007 data 
from |H10| Although our lightcurve is much more precise, we did 
not record enough periods to improve the result from ,H10t 

The lightcurve exhibits a clear asymmetry. This asymme- 
try can actually also be guessed by eye in the HIO" data where 
the oscillation maximum seems slightly off'set from a purely si- 
nusoidal signal (Eig. 6 of their paper). In order to confirm this, 
we have processed the publicly available SuperWASP WASP- 19 
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Fig. 2. (Top) 23-days lightcurve of the WASP- 19 star showing its variability (shifted for clarity). The stellar variability is over-plotted in red, 
including a final filtering step to remove residual systematics (sliding median of 6.25 hours in width). (Bottom) Same as the top lightcurve, but 
with the stellar variability removed (see text for details). 



data and performed several tests showing that the asymmetry is 
present in 2007 and comparable in shape (or trend) to what we 
obtained in 2010. In 2008, we find a change in the asymmetry, 
although with an amplitude about twice fainter than in 2007 and 
2010. This appears to result directly from the presence of stellar 
spots and their progressive evolution. 

3.2. Primary transit 
3.2.1 . Transit parameters 



Table 2. Best fit transit parameters from the ASTEP 400 observations. 



The phase-folded transit on Fig. |4] was fitted using the |MandeI 
et al.| ( |2002| ) analytical transit formulation and an IDL imple- 
mentation of the Levenberg-Marquardt algorithm (LMFIT) to 
find the best fit (black curve on Fig.|4]). We assume a circular 
orbit to define the timing of the secondary eclipse in line with 
the previous secondary eclipse measurements. The stellar limb 
darkening was adjusted with a quadratic law. The obtained tran- 
sit parameters of WASP- 19b are summarized in Table [2j where 
P is the orbital period, Tq the reference mid transit time, Rp the 
radius of the planet, 7?^ the radius of the star, a the orbit semi- 
major axis, / the orbital plane inclination, b the impact parameter, 
7i and 72 the quadratic limb-darkening coefiftcients, e the eccen- 
tricity, and oj the argument of periapsis. 71 and 72 were derived 
from previous studies ( Hebb et al.|2010[ Hellier et al.|2011|| ) and 
tabulated data from [Claret & Bloemen| ( |2011| ). However, the de- 



Parameter 



Value 



Error 



P[d] 
To [HJD] 

a/R^ 

in 

b 

n 

72 

e 



0.788840 
2455317.26998 
0.141 
3.560 
79.6 
0.642 
0.42 
0.29 





(fixed) 
±0.00006 
±0.001 
±0.036 

±0.3 
±0.023 
(fixed) 
(fixed) 
(fixed) 
(fixed) 



pendency of the fit on these parameters is very low at this level 
of photometric precision. 

3.2.2. Fluctuations in transit depth and timing 

The dispersions in transit length, timing, and depth, are compati- 
ble with the level of pink noise, assuming no intrinsic variations. 
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Fig. 3. Summary of noise statistics for each night computed on the lightcurve residuals (data subtracted from the model) for each complete run 
(or equivalently for each observing night). Labels 'P', 'R' and 'W stand for pink, red and white, respectively. The run (or night) number as well 
as the number of considered data points {np) is indicated above each frame. The abscissa correspond to HJD-2455000. 



The primary transits were first fitted altogether to define 



mean values (see Fig. |B.l| in Appendix|B]). We then allowed the 
depth, timing or duration to vary between the individual transits. 
Figure|5] shows the resulting dispersion between transits. The 
uncertainties are calculated from the pink-noise of the respective 
nights derived above. 

We find that the fluctuations in timing and duration are com- 
patible with the error estimates, indicating that our treatment of 
correlated noise is satisfactory. The point-to-point dispersion of 
the depth is also compatible with the error estimates. 

3.3. Secondary eclipse and phase effect 

The secondary eclipse is detected on the folded lightcurve, 
where an approximately synchronous phase efl'ect is also evi- 
denced (inset of Fig.|4]). In order to determine whether to con- 
sider all of the observing night, or only specific good ones, we 
analyzed the behavior of the noise at the eclipse location for an 
increasing number of N nights. In order to determine the good- 
ness of nights, we introduce the eclipse noise cr^^. 



t 



\ 



10 15 
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Fig. 5. Fluctuations of the phase timing (left panel), length duration 
(middle panel) and depth (right panel) of individual transits of WASP- 
19b relative to the values from the combined transits. The crosses cor- 
respond to runs (nights) with no transit event, and orange point to nights 
with incomplete primary transit. 



+ o-z 



where the (Tw^ and cr^^ values are taken from 
mated for each observing night at the eclipse location, thus tak- 



Fig.[3] 



ing into account the number rik of data points at the bottom of the 
eclipse event. The cumulative noise cr^^ in the folded lightcurve 
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Fig. 4. Phase-folded lightcurve of WAPS-19b, where the gray and blue dots correspond to the original data, and to the binned data at 0.001 phase 
units (blue), respectively. The black line is the best fit to the transit (see Table [2] for detailed primary transit parameters). The inset frame is the 
normalized flux scale zoom on the secondary eclipse (full phase coverage in abscissa) with a bin size of 0.01 in phase - see Sect. |3.3| for detailed 
secondary eclipse parameters. 



is then defined by, 



1 



2 2 



with 



(Tej^ is plotted on Fig.|6j and shows that the eclipse noise always 
improves with increasing number of nights. The plot also dis- 
plays the cumulative white noise, as well as the value of cr^^ for 
each additional night (in a non-cumulative way). 

The lightcurve shown in fig. |4] was obtained from all the 23 
nights, but we checked that the actual measured features did not 
originate from the few, noisier observing runs only. 

After pre-processing the data, we obtain a secondary eclipse 
fit with a 390 ± lOOppm depth, using the same Levenberg- 
Marquardt algorithm as in Sect.|3.2.l] Note that the ilOOppm 



error only corresponds to the model fitting uncertainty, and 
should not be confused with the uncertainty on the real depth 
of the eclipse (detailed below). Considering previous works at 
other wavelengths ( Gibson et al.| ( |2010| )) we performed the fit of 
the secondary with a fixed phase of 0.5 consistent with a circular 
orbit. 

In order to probe the reality of this detection, we performed 
several tests. Based on the measured points themselves, we car- 
ried out a Monte-Carlo simulation consisting in adding artificial 
secondary eclipse events to the actual data. These events have 
a random depths as well as random periods (excluding a range 
containing the primary transit period). Random periods rather 
than random phase with the WASP- 19b period was used in order 
to scramble the fluctuations of the noise (especially the red com- 
ponent) as well as systematics. For each of these new data sets, 




1200 




Fig. 6. Standard deviation of the folded lightcurve for 1 hour bins 
against the number of folded nights (from best to worst). The numbers 
along the curve correspond to each of the newly added night number: 
the first point corresponds to a single night whilst the latest corresponds 
to the 23 nights. The lower panel shows the cTg^ value for each night. 
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Fig. 7. (top) Result of the Monte-Carlo simulation points (see text for 
details). The blue circles show the detected events (i.e. with the phase 
of the fit matching the expected eclipse location +2.5% in phase unit), 
whilst the red dots show non detected events (i.e. the result of the fit not 
matching the expected phase location +2.5%). (bottom) The detection 
(blue) and false alarm (red) probability curves derived from the Monte- 
Carlo simulation. The detection probability - for a given depth - is the 
ratio between the number of detected events (blue circles) over the total 
number of events (red plus blue). 



an eclipse detection fit is performed so that we can determine 
the number of detected events, i.e. the detected events which 
phase matches the introduced one, without constraints on the 
depth (what we are looking for). The tolerance for the detected 
phase is ±2.5%, implying a false event detection probability of 
5% assuming that these random false events are uniformly dis- 
tributed in phase. This tolerance (±2.5%) was arbitrarily set as 
the 2.5-cr dispersion of the fitted phases at an eclipse depth of 
500 ppm. The results from this simulation are shown on Fig. [7] 
where 20000 realizations have been drawn. This allows us to 
assess: 

• what we call the detection sensitivity to known, injected 
events of a given depth, which is computed as the ratio of the 
number of detected events over the total number of injected 
events at this given depth (i.e. a statistic out of an horizontal 
cut of the top panel of Fig.|7]), 

• and the ratio of detected (or spurious) events for a given mea- 
sured depth, which is computed as the ratio of the number of 
detected events over the total number of measured events at 
this depth (i.e. a statistic out of a vertical cut of the top panel 
ofFig.|7]). 

For injected eclipses with a depth 390 ppm, the detection sensi- 
tivity at this depth is 69% (with a bin width of 40 ppm), while 
the probability of detecting an eclipse event at this depth is 43%. 
However this latter number does not take into account the prob- 
ability that these events are due to false events. 

The uncertainty on the real depth of this event is simply de- 
termined by building an histogram of the injected transit depths 
at the measured depth of 390 ppm. The result is shown on Fig.[8] 
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Fig. 8. Histogram showing the depth of injected events that were 
detected in the MC simulation, at the specific measured depth of 390 
(a vertical +40 ppm bin through the diagram of Fig.[7]-top). The blue 
curve shows the best-fit Gaussian of the histogram. 



In order to quantify the average depth and error, we choose to fit 
a Gaussian to the obtained histogram. The best fit curve gives 
a depth of 391 ± 5 ppm and a cr = 194 ± 4 ppm, which we ac- 
count for the depth and true uncertainty of the eclipse. We adopt 
the rounded values of 390 ± 190 ppm as the final result for this 
work. For the phase of the secondary eclipse, a similar process 
was applied and lead to a shift in phase (relative to 0.5) found 
by the model fitting of 45 ± 15 s, while the error estimate on the 
"real" phase oflfset is around ±30 minutes. At this level of preci- 
sion, our result is matching previo us works ( [Gibson et al.|2010| 
Hellier et al.|2011[ [Anderson et al.|2011) that find no significant 
phase shift, or one well below our estimation precision. 

The inset in fig. |4] exhibits a clear phase-eflTect that is fitted 
with a 308 ± 90 ppm peak-to-valley amplitude for a sine func- 
tion, and with a phase ofifset of -0.04 ± 0.03 relative to the 0.5 
phase. The lightcurve was obtained by filtering the original one 

The 



to get rid of the stellar variability as described in Sect 2.3 



fact is that the last filtering step would introduce an attenuation 
of the phase eflTect at this period (if any) by a factor that we esti- 
mated to 3.2 ± 1.2 using a method similar to that of the secondary 
eclipse Monte-Carlo test. In order to assess the possible origin of 
this phased function, and understand whether it is real or caused 
by trends and/or the filtering process, we performed two tests: 
trying random combination of 10 nights to generate the folded 
lightcurve, and worked on synthetic data. In the first case, the 
phase function appeared almost in all cases, so that we are confi- 
dent that the eflTect is not caused by specific nights. In the second 
test, we injected random trends and phased signals to a synthetic 
lightcurve (same time sampling as the original data) so as to 
mimic the measured folded lightcurve. These tests suggest that 
this phased oscillation is compatible with a mix between residual 
from strongly attenuated trends and an attenuated phased signal 
that appears to be real (be it of astrophysical origin or caused by 
systematics), that cannot be clearly disentangled. However, we 
find that if there is a phased signal, its amplitude would have to 
be larger (> 500 ppm) than the planet phase eflTect alone to cause 
such a well-in-phase result. In conclusion, and in the absence of 
other work trying to measure the phase eflTect, we can only be 
very cautious about its reality in our data, but the correct phas- 
ing and symmetry of this hypothetical phase eflTect is somehow a 
bit disconcerting. Our observation look promising regarding the 
possibilities to assess such kind of eflTect in upcoming observ- 
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ing campaigns. Note that we also fitted the lightcurve letting the 
phase eff'ect and eclipse phases as free parameters, as a test of 
robustness of our secondary eclipse and phase detection. The 
best fit gives a secondary eclipse depth and phase of 374 ppm 
and -0.007, and a phase eff'ect amplitude and phase of 300 ppm 
and +0.05. 



4. Analysis 

4.1. Brightness temperatures and albedo 

The analysis of the spectra of exoplanets is complicated by their 
unknown chemical composition, cloud structure, their inhomo- 
geneous temperature field and complex dynamics. The detec- 
tion of a secondary eclipse is a measurement of the emission of 
the day side of the planet, including both thermal emission and 
reflected light. Disentangling the two generally requires an ex- 
tremely accurate full phase lightcurve at visible wavelengths, so 
that both the brightness temperature of the day side (during sec- 
ondary transit) and the temperature of the night side (during pri- 
mary transit) can be measured, as exemplified for TrES-2b and 
HAT-P-7b in the K epler field ( Kipping & Bakos] [2QTTt |Chris- 



Table 3. Measured secondary eclipse depths and brightness tempera- 
tures from eq. ([TJ as a function of wavelength for WASP- 19b. 



|tiansen et al.|2010| ). Unfortunately, the ASTEP phase lightcurve 
of WASP- 19b does not reach the accuracy required to break the 
degeneracy. The analysis must therefore rely solely on estimates 
of the day side brightness temperatures at various wavelengths 
and comparison with other cases. 

The day side brightness temperature Tbright of the planet is 
calculated from the depth d of the secondary transit following 
Cowan &AgoI| ( [20TT] ): 



d = 



(1) 



where Rp and R^, are the planetary and stellar radius, respec- 
tively, 7r^(rbright) corresponds to the assumed blackbody emis- 
sion flux of the planet's day side, and FiT^f^^^) is the stellar flux. 
In this work, we approximate the latter by a theoretical stellar 
spectrum of the same efl ^ective temperature and metallicity from 
IHauschildt etaL] ( |T998l ). 

The secondary eclipse depths measured at difl'erent wave- 
lengths and the corresponding brightness temperatures calcu- 
lated with eq. ^ and for T^^ = 5500K, R, = 0.93 R©, R^ = 



1.3 Rjup (HIO) are presented in Table B.l Clearly, the measure- 
ments at short wavelengths correspond to smaller depths and are 
thus more challenging. The corresponding brightness temper- 
atures are always very high and indicate that the day side of 
WASP- 19b is very hot. For comparison, the planet's equilib- 
rium temperature for an assumed Bond albedo (assuming the 
incoming stellar flux is redistributed uniformly on the planet's 
surface - see Sau mon et al.|1996 ) is Tq^^ = T^f^^^^DIAR^yi^ = 
2000 ± 40 K, where D is the star to planet distance. The maxi- 
mum efl'ective temperature that would be emitted by the day side 
if it was a blackbody with no redistribution of heat from the day 
side to the night side is rinaxirr,o = (8/3)^/^req,o = 2556 ± 51 K. 
Brightness temperatures in excess of that value are possible at 
selected wavelengths because of opacity variations or in the vis- 
ible because of reflected light. In the present case, it is interest- 
ing to note that the measured temperatures are slightly in excess 
but still compatible with that maximum efl'ective temperature in 
bands that span a wide wavelength range, from 0.42 to 2.2 jam. 

he geometrical albedo that would be necessary to explain the 
measurements in the visible isAg = d (D/R^)^, i.e. in the case of 
WASP-19b in the ASTEP 400 band, we obtain = 0.27 + 0.13. 



Wvl. range (|Ltm) 


Depth (ppm) 


T'bright (K) 


References 


0.42-0.95 
0.81-1.01 
1.525-1.775 
1.965-2.215 


390 + 190 
880 + 190 
2760 + 440 
3660 + 670 


2690^1^0 
2640+^^ 
2740!|» 
2680![^o 


1 
2 
3 
4,5 


3.225-3.975 


4830 + 250 


2330!f, 


5 


4.0-5.0 


5720 + 300 




5 


5.1-6.5 
6.55-9.45 


6500+ 1100 
7300 + 1200 


2244+230 
2240!2| 


5 
5 



References. (1) This work (see Fig.f T]); (2) [Burton et al.| ([20T2l); (3) 
Anderson et al., ( .2010j ; (4) Gibson etaTlpOlOl ); (5) [Anderson et al.| 
l20Tl1 ). 



This is to be regarded as an upper limit on the geometrical albedo 
at these wavelengths. Because of the similarity of the values 
of the brightness temperatures in the visible and near-infrared, 
it seems however more likely that the observed signal is due to 
thermal emission from a very hot day side, with a very ineflftcient 
redistribution of heat from the day side to the night side. 



In order to put the WASP- 19 observations in context, we 
plotted in Fig. [9] the exoplanets which show a significant (> Icr) 
detection of a secondary eclipse both in the visi ble and in the in - 
frared, as taken from the compilation of, Cowan & Agol[ ( [2011 ). 
We use the models of [Fortney et al.[ ( [2005 [2008| ) tocompare 
to the observations of these planets. These are strictly one- 
dimensional day side models that allow a certainly fraction of 
absorbed light to be "lost" to the night side, which allows the 
day side Tqq to be adjusted. Here the model fluxes are converted 
to brightness temperatures. Except for CoRoT-lb, all models in- 
clude TiO and VO as day-side absorbers, which yields day-side 
temperature inversions. The high K-band fluxes for CoRoT-lb 
suggest a high temperature in this water opacity window, which 
argues against a temperature inversion. Solar abundances and 
equilibrium chemistry are used for all models and no ad-hoc 
adjustments to chemical mixing ratios were made to yield im- 
proved fits. For the models shown in Fig.|9] the red curves are 
due to thermal emission, and the green curves are thermal emis- 
sion plus reflected light. The green curves turn dramatically up- 
ward at short wavelengths due to an increase in the reflected light 
component, but the geometric albedos at these wavelengths are 
still quite low. 



Clearly, for a number of the planets shown here, the sim- 
ple models do not capture the true complexity of these atmo- 
spheres. The models generally have redistribution over the day 
side only, leaving only -0-20% of the absorbed energy for the 
night side. For WASP- 19b, the model has a day side reff=2300 
K. The planet's true day side Tqq must however be larger since 
the model is not bright enough from -0.8-3 jim. In the model 
we find the optical emission from the planet with ASTEP 400 is 
dominated by thermal emission. 
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in the star. As shown for the OGLE-TR-56 system by Carone & 
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Fig. 9. Ratio of the brightness temperature to the maximum irradia- 
tion temperature (at the substellar point - see text) versus wavelengths 
for exoplanets with detected visible and infrared secondary eclipses. 
Green and red curves correspond to the total modeled spectrum and to 
the thermal emission only, respectively. The ellipses account both for 
the bandpasses of the observations (horizontally) and the uncertainties 
of the measurements (vertically). The colors indicate whether the data 
originate from space telescopes (blue), mid-latitude ground based facil- 
ities (orange), or ASTEP 400 (red). 



4.2. Tidal evolution 



As shown in|H10', WASP- 19b is at risk of engulfment due to its 
extreme proximity from its central star. Specifically, the inward 
migration timescale for a planet of mass Mp in a circular orbit of 
period P around a star of mass M*, radius and spin period P^, 
is Parker & Ogilvie|2009| : 



12.0 Ma 



9±. 
106 



Mo 



8/3 



X 



M 



Jup 



-1 



Ro 



-5 



1 day 



P 



-1 



where Q'^ is the equivalent tidal dissipation factor in the star (de- 
fined as 3/2g^/^2,*, where 2* is the standard tidal dissipation 
factor and ^2,1 its Love number - see Barker & Qgilvie||2Q09| ). 
Using the parameters from Table|2]we obtain: 



3.8 Ma 



106/' 



(2) 



This implies that either we are extremely lucky to see a system 
in its few last million years of existence, or more likely the tidal 
2^ is significantly larger than the assumed 10^ typically required 



to explain the circularization of binary systems (e.g. [Meibom & 
[Mathieu 2005). 

The spin period of the star derived from both HIO and this 
work enables to put constraint on the amount of tidal dissipation 



Patzold (2007), the tides raised by the companion on the cen- 



tral star yield a competition between the orbital migration of the 
companion, the magnetic braking of the star, and its possible 
tidal spin-up. As a result, values of 2^ that are too low may be 
ruled out because they lead to a spin-up of the star that is too 
efficient. For the OGLE-TR-56 system, Carone & Patzold find 
that 2; > 3 X 10^. 

We apply the same method to the WASP- 19 system, using 
planar -dynamic al evolution equations for the system (B arker^&| 
Qgilvie 2009; Guillot & Havel||2QTT] ), stellar evolution models 



([Morel & Lebreton|2008 ) anTa magnetic braking relation from 
We explore the solutions at constant 2* 
We assume an arbitrary initial eccentric- 



Bouvier et al, 



between 10" andW, 
ity of 0.7, but even at the highest values of 2* considered, we 
find that the solutions are independent of this assumption due 
to the very rapid circularization of the orbit in the first few mil- 
lion years of evolution. The stellar spin is initialized at 4 days , 
a value representative of T-Tauri stars (see |Bouvier et al.|1997 ). 
The solutions are found to be weakly dependent of this param- 
eter. For each assumed value of 2*. we try to find the value 
of the initial semi-major axis of the planet that yields a solu- 
tion in agreement with our observations. The results are pre- 
sented in Figure[To| Only values of 2* larger than about 3x10^ 
yield both an orbital period and a stellar spin period that match 
the observations after about 1 to l.SGa of evolution. The lat- 
ter ages correspond to the lower values of 2* and effectively 
require a spin-up of the parent star, as indicated by the compari- 
son with other exoplanetary systems and stellar isochrones (see 



Husnoo et al.|2012|). These results are in agreement with the val- 
ues obtained in |H 10 with slightly different assumptions. While 
a more extended exploration of the parameter space would be 
desirable but beyond the scope of this article, we therefore feel 
that our results are relatively robust. Surprisingly perhaps, the 
value obtained is very similar to that obtained for OGLE-TR-56 
by Carone & Patzold. While WASP- 19 is a G8V type star with 
a mass slightly smaller than our Sun, OGLE-TR-56 is an F-type 



star, with an effective temperature around 6100 K (Torres, et al. 



|2008| ). In both cases, the companions have masses slightly above 
Jupiter's (1.35 Mj^p for the latter), indicating the dissipation pro- 
cess is weak for Jupiter-mass planets, in line with theoretical 
predictions of dissipation by inter nal gravity waves ( [Ogilvie & 
Lin|2007t [Barker & Ogilvie|2010| ). 



5. Conclusion 

We have presented observations of WASP- 19 with ASTEP 400 
obtained during a 24 days observing period through May 2010 
that show that extremely stable and precise visible photome- 
try and near-continuous observations are achievable from the 
Concordia station, Antarctica. The data appear to be of excel- 
lent quality, with hour-timescale pink noise averaging 0.7 mmag, 
comparable to very good ground-based transit photometry. The 
observations confirm that WASP- 19 is a relatively spotted star, 
with ~ 1 % peak-to-peak photometric variations for this period of 
observation, and a spin period of 10.5 ± 0.2 days. 

A detailed analysis of the transit signal yields a detection of 
a dip of 390 ± 190ppm at phase 0.5 consistent with the pres- 
ence of a secondary eclipse. It would correspond to a brightness 
temperature of the day side of the planet of 2690^220 ^^^^ 
slightly in excess of the maximum possible effective temperature 
for no redistribution of heat of 2556 ± 51 K, but well in line with 
brightness temperatures obtained in the and K bands. This 
seems to imply that the emission is due to a very bright day side 
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Fig. 10. Dynamical evolution of the WASP-19 star-planet system 
showing the orbital period versus the stellar spin period, in five diff'er- 
ent cases assuming values of the stellar tidal dissipation factor Q'^ (see 
text) ranging from 10^ (2 cases) to 10^. The corresponding ages, in mil- 
lion years, are shown by the color scale to the right. The observational 
constraints are represented by a square. Only models with values of 
above 3. 10^ are able to match the observations. 



and that heat redistribution is very inefficient. Alternatively, the 
signal observed in the visible may be due to direct reflection of 
incoming starlight, in which case we can derive a maximum geo- 
metrical albedo Ag = 0.27 + 0.13. Unfortunately, no atmospheric 
model appears to correctly account for all observations of the 
planet so far. 

Separately, the ability to fully characterize the system, and 
in particular the stellar spin period, allows us to put useful con- 
straints on the tidal interactions between the star and its planet 
and specifically on the tidal dissipation in the star Q'^ > 3 x 10^, 
a value much larger than that required to explain the circulariza- 
tion of binary stars. The detailed study of close-in giant planets is 
thus extremely fruitful both to understand planetary atmospheres 
and tidal interactions between stars and their companions. The 
ASTEP 400 detection of a secondary transit in the visible is the 
first achieved from the ground. We also identified a phase func- 
tion that may originate from the combination of residual trends 
and a real periodic signal that cannot be clearly disentangled 
at that point. Further investigations could be carried out with 
ASTEP 400, but would require observations for a significantly 
longer duration. 



Appendix A: OlS photometric correction 

Appendix A.1: OlS principle 

We recall the OIS principle that consists in subtracting a refer- 
ence image (R) to the current science one (/). This reference 
image is the best one in the time series, in the sense of sharpness 
and sky background level. It is then convolved to an optimal 
kernel which depends globally on the current science image, and 
locally to the position in the image (space- varying kernel). Once 
the reference image is convolved by the optimal kernel (K), it has 
the same average intensity and PSF shapes as the current science 
image. After subtracting the convolved reference and the sci- 
ence image, we obtain an image with PSF subtraction residuals, 
whose flux need to be measured by aperture photometry. This 
operation can be written as. 



where Fref is the flux computed on the reference image with the 
aperture A. With our notation, we could write F^^f = AR. 

Appendix A.2: Choice of the aperture radius 

The radii of the apertures A (which are diflTerent for each star 
in the reference image) depend on the intensity of the stars and 
on the global FWHM of all PSFs. The aperture is allowed to 
vary for each star so as to minimize the global photometric noise 
structures in the lightcurves. This is done at the expense of a 
photometric bias - as explained below - that can be compensated 
using a statistical approach detailed in the next section. The net 
result is a significant gain in SNR. 

In order to optimize the photometric SNR, a simple method 
consists in selecting an aperture radius by determining all pix- 
els with a level higher than (for example) three times the back- 
ground noise. By this process, each star in the reference image 
has its own aperture radius A. Note that fainter stars will have 
smaller reference aperture radii by using this method. 

These apertures will give us the reference fluxes F^Qf for each 
star, but if we use these apertures directly on the subtracted im- 
age, they will not take into account the PSF FWHM broadening 
(it is necessarily a broadening since the reference image is cho- 
sen as the sharpest image). A way to account for these FWHM 
width fluctuations is to convolve the circular reference aperture 
by the local kernel K. The flux is then expressed as. 



F = Fref + 



(A <^K)[I-(R<S> K)] 



(A.2) 



Apart from broadening the apertures proportionally to the 
current science FWHM, this aperture convolution has the bene- 
fit of 'apodizing' the resulting aperture, thus reducing the back- 
ground noise contribution. To better understand this eflTect, let 
us consider the case of a star whose photometry is perfectly 
constant over time. Its reference flux is determined on a ref- 
erence image with a given circular aperture radius. In another 
image where the PSF FWHM is larger, the star's flux will re- 
main the same, but will be spread-out because of degraded see- 
ing. If we use the circular aperture convolved with the kernel K 
as the new aperture, we can understand why the outer pixel val- 
ues (that should contribute to the true photometry) will be artifi- 
cially decreased, whereas the background noise, in areas where 
it becomes dominant compared to the star signal, will be further 
attenuated. 

Thus the image subtraction algorithms intrinsically and glob- 
ally produces more accurate (cleaner) lightcurves than aperture 
photometry, but has a tendency to underestimate the photome- 
try (and the noise level) in the case where the above-described 
aperture convolution procedure is used. 

It should be pointed out that the photometric bias is more 
important for fainter stars due to the combination of two eflTects: 
firstly, as mentioned above, the associated reference image aper- 
ture radius will be smaller than for a brighter star, and secondly 
because the kernel size is approximately constant over the entire 
image. Since the average kernel size directly depends on the see- 
ing, the photometric underestimation bias will be even stronger 
for fainter stars during bad seeing episodes, and can reach up to 
about 30% in our data. 

In order to compensate for this photometric bias, we need to 
introduce a correction factor a such that. 



F = Fref + 



A [/-(/? K)] 



(A.1) 



F = Fref + «f 



(A (S>K)[I-(R^ K)] 



(A.3) 
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Fig. A.l. Plot of the reference image fluxes computed with a circu- 
lar aperture (Fref) versus the same reference image computed with the 
kernel-convolved aperture A <S) K (Fcref). The fitted curve is shown in 
red, as well as the ideal linear trend (dashed line). 



Appendix A.3: Photometric correction factor a 

To assess the a factor, we simply plot the reference fluxes Fref of 
each star against the reference PSFs, convolved by the modified 
aperture, that is Fcref = (A^ K){R^ K)/ \\ K f. This shows 
how the reference fluxes (in the reference image) compare to the 
fluxes when we use the convolved aperture A <^ K. If the field 
contains enough stars, a fitting curve can be adjusted in order to 
give the photometric compensation a for any measured photo- 
metric residual (Fig. |A.l| ). The photometric correction factor a 
needs to be computed for every image since the optimal kernels 
are changing with each image. The uncertainty on a was evalu- 
ated by scrambling all the points, randomly switching them (or 
not) with its immediate neighbor, then performing the fit and re- 
peating the procedure 1000 times. This test shows that a is stable 
at the level of about 1 % up to magnitude R = 16, and raising to 
about 5 % at 7? = 18. An additional test showing the robustness 
of the method is that the correction factor a is very well corre- 
lated with the PSF FWHM (thus consistent with the explanation 
given above). Although we cannot verify the exactness of the a 
factor for each of the stars, it appears to be an efficient first order 
photometric correction. Most importantly it does not introduce 
spurious noise structures in the lightcurves as it is the case of 
WASP- 19 (see below). 



a small crop (600 x 600 pixels) rather than on the full image to 
increase the OIS precision (recalling that the kernel is allowed 
to vary according to second order variations). The WASP- 19 
field is not as crowded as other of our fields where the correc- 
tion relies on much larger star counts. We therefore computed 
the photometric compensation factor for each image individu- 
ally, and used a temporally smoothed version of it to compensate 
for the original OIS data. In that way, the original signal to noise 
is preserved, while correctly compensating for the photometric 
bias, mostly visible when deriving the primary transit parame- 
ters. Note that without compensating for the 10% photometric 
bias, the obtained primary transit parameters are way beyond 
the ones obtained with the compensated data (including error 
bars). The compensated data produce parameters that are totally 
consistent with v alues derived i n previous papers (e. g. [Ander^ 
son et al.| ( [2QTQl )). Especially in'Hell ier et"aL] pOTT ) a WASP- 
19b transit lightcurve obtained with the NTT in the Gunn r-band 
is shown (see their Fig. 1 -bottom), that was recorded just three 
months prior to our own observations, and is totally consistent 
with our own results. 



Appendix A.5: Comparison between aperture photometry 
and OIS 

Figure [Al2| is a comparison between the residuals obtained for 
both the aperture photometry and the OIS algorithms in the 
whole WASP- 19 field. The lightcurve calibration procedure is 
identical for both algorithms (only the raw lightcurves are diff'er- 
ent). 



Appendix AA: The case of the WASP-IQ field 

In the case of WASP- 19, the transit photometry was underesti- 
mated by approximately 10%. We also found that the described 
photometric correction procedure can be the source of additional 
noise (specifically for that field) if the correction factor is com- 
puted and applied to each individual image. Indeed, in that case, 
the correction is not simply a scaling factor, but seems to gener- 
ate more high frequency fluctuations (typically of about an hour 
length). The reason is that for WASP- 19, we applied the OIS on 



Appendix B: Day-to-day transit lightcurves 
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Table B.l. Basic data and primary transit parameters for each of the observing night. 



Night 


(} data 




hA 




Pinknoise 


Whitenoise 


-K-^dnoise 






(HJD) 


(min) 




(ppm) 


(ppm) 


(ppm) 


2 


257 


2455318.05715 


106.25 


0.150 


722 


429 


580 


3 


275 








434 


392 


186 


4 


273 


2455320.42864 


76.44 


0.135 


661 


439 


494 


5 


284 




99.19 


0.145 


594 


342 


486 


6 


293 


r\ A c C ^ r\r\ f\r\r\ a 

2455322.00248 


95.50 


0.139 


484 


429 


224 


7 


24 














8 


217 








769 


426 


640 


9 


286 


2455325.15851 


100.06 


0.135 


776 


374 


680 


10 


230 


2455325.94792 


94.95 


0.134 


683 


571 


375 


11 


286 








415 


342 


236 


12 


341 


2455328.31349 


96.78 


0.139 


734 


326 


657 


13 


311 


2455329.10318 


95.34 


0.140 


962 


409 


871 


14 


260 








730 


370 


629 


15 


351 


2455331.46889 


108.73 


0.155 


728 


349 


639 


16 


295 


2455332.25741 


91. lA 


0.140 


646 


332 


554 


17 


345 


2455333.04710 


106.19 


0.147 


595 


340 


489 


18 


284 


2455333.83724 


99.80 


0.144 


708 


412 


576 


19 


360 


2455318.05715 


106.25 


0.150 


441 


344 


275 


20 


127 








741 


363 


646 


21 


367 


2455336.99102 


95.57 


0.139 


591 


372 


459 


22 


359 








499 


338 


367 


23 


377 


2455339.35751 


95.94 


0.139 


955 


389 


872 


24 


338 


2455340.14610 


95.77 


0.133 


647 


435 


480 
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Fig. B.l. Day-to-day lightcurves of the primary transits events, along with the best-fit model (from parameters of Table[2]) overplotted in red, and 
the residual points, shifted in ordinate for clarity. 
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